Dynamic equilibrium Mechanism for Surface Nanobubble Stabilization 
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Recent experiments have convincingly demonstrated the existence of surface nanobubbles on 
submerged hydrophobic surfaces. However, classical theory dictates that small gaseous bubbles 
quickly dissolve because their large Laplace pressure causes a diffusive outflux of gas. Here we 
suggest that the bubbles are stabilized by a continuous influx of gas near the contact line, due to 
the gas attraction towards hydrophobic walls (Dammer & Lohse, PRL96, 206101 (2006); Zhang et 
al, PRL98, 136101 (2007); Mezger et at, J. Chem. Phys. 128, 244705 (2008)). This influx balances 
the outflux and allows for a meta-stable equilibrium, which however vanishes in thermodynamic 
equilibrium. Our theory predicts the equilibrium radius of the surface nanobubbles, as well as the 
threshold for surface nanobubble formation as a function of hydrophobicity and gas concentration. 



Various recent studies have revealed the existence 
of nanoscopic soft domains at the liquid-solid interface 

0, i, i i a i a a a 0, [n e 0, q m, \m- 

As atomic force microscopy (AFM) shows, these soft do- 
mains resemble spherical caps with heights of w 10 nm 
and diameters of 2R « 100 nm, corresponding to a small 
contact angle (on the gas side) of 9 ^ 10°, see figure 
1 for a sketch. The most consistent interpretation [l3| 
of these soft domains is that they are surface nanobub- 
bles, i.e., nanoscale gas bubbles located at the Hquid-sohd 
interface. Spectroscopic studies ^^'^ X-ray reflec- 

tivity measurements [Tol . [20l [2l| experimentally confirm 
the presence of gas at the wall. Moreover, the size and 
density of these objects depend on the dissolved gas con- 
centration and they disappear completely when the liq- 
uid is degassed [11 [H [B, [H HI] . CO2 bubbles have 
a much shorter lifetime (only 1-2 hours) as compared to 
air bubbles, due to the much lower pressure of CO2 in 
the atmosphere and its better solubility in water p^ . 
The existence of these surface nanobubbles has poten- 
tially great technological importance, as they have been 
shown to result in fluid slip in hydrophobic surfaces, lead- 
ing to a potentially lar ge r eduction in fluid dissipation for 
flows in small devices [2I, [H, . 

Observations of long living surface nanobubbles fl atly 
contradicts the classical theory of bubble stability [2a |. 
Small bubbles have large internal gas pressures, in or- 
der to balance the compressive action of surface tension. 
For a bubble in water of radius R — 50nm, which with 
above contact angle 9 = 10° corresponds to a radius of 
curvature Rc = R/ sin0 ss 250nm, and surface tension 
7 — 73mN/m (20° temperature), the gas pressure ac- 
cording to classical macroscopic theory (neglecting the 
disjoining pressure [131) is Pgas = "^1 / Rc = 2jsm9/R = 
0.58MPa. In employing Laplace's equation we also ne- 
glect the bubble-stabilizing contribution from hydroxide 
ion adsorption on the bubble surface [2^] . If the resulting 
charge density is constant, this ion adsorption contributes 
a radius- independent term to Pgas, which is ncgiibible 
for small enough bubbles. Henry's law then dictates 
that the concentration of gas at the bubble surface is 



c(i?) = coPgas/Po, where cq is the saturated concentra- 
tion of gas at atmospheric pressure Pq. This is necessarily 
much larger than the gas concentration Coo < cq far away, 
leading_^to a diffusive outflux of gas and bubble dissolu- 
tion 13 ■ V the contact angle 9 remains finite for 
R ^ and assuming that Laplace's law remains valid at 
molecular scales, then c{R) formally diverges, signalling 
the breakdown of the continuum approach. 

Various suggestions have been made to explain the long 
lifetime of surface nanobubbles. Among these, the reduc- 
tion of surface tension for large curvature on small scales 
[sol . [3l| (operative for bubbles smaller than ^ 5nm) , the 
oversaturation of the liquid around the nanobubbles with 
gas [3, [l3|, the effect of induced charges in the Debye 
layer around the bubble interface [12], or the stabiliza- 
tion of bubbles due to contaminant molecules in the water 
[33| . similarly as has been suggested for microbubbles in 
contaminated water [34| . In this paper we propose an 
alternative mechanism for bubble stabilization, namely, 
that the gas outflux is compensated by gas influx at the 
contact line. 

The diffusive mass outflux from a surface nanobubble 
is given by 

Jout{R)^Dj dAn-\/cKi2TTD rdrd^c, (1) 

where D is the diffusion constant of gas in the liquid 
(typically 10~^m^/s) and n is the normal vector to the 
droplet surface. We have assumed that the bubble is 
nearly flat, so that the diffusive gradient is primarily in 
the z direction, perpendicular to the solid surface. This 
diffusive flux can then be evaluated by solving the steady 
state diffusion equation in the liquid. The result is 

JoutiR) = TiRD{c{R)-Coo) - ^RD _ Co 

(2) 

The volume flux rate jout[R) — Jout{R)/c{R) is then 
given by 



j,^,iR) = :^^^^nRD(l- ' 



c{R) 



c{R). 



(3) 



The volume of the bubble thus decreases linearly with 
time, typically with a volume flux rate sa nRD « tt • 50 • 
10~^ • 10~^m^/s, when assuming complete degasing. For 
an initial bubble of radius 50nm, this implies a dissolution 
timescale of ~ 1/iS. 

Stabilization against dissolution requires a physical ef- 
fect to cancel this diffusive outflux. Although intermolec- 
ular forces might modify the liquid surface tension from 
its macroscopic value [l^, HH , the energetic cost of cre- 
ating surface energy ensures that the gas pressure in the 
bubble is always higher than that of the surrounding liq- 
uid, and hence diffusive outflux necessarily persists even 
when interactions with the solid surface are accounted 
for. Stabilization can be achieved however, by mecha- 
nisms causing an influx of gas into the bubble. Such an 
influx need not occur uniformly across the bubble sur- 
face, but can be spatially concentrated. In particular, it 
is quite natural to consider influx mechanisms near the 
contact line, where intermolecular forces are most signif- 
icant. For a sketch of the gas flow directions, we refer to 
figure 1. The magnitude of any contact line dominated 
influx increases with the circumference of the bubble, and 
hence could compete effectively with the diffusion medi- 
ated outflux eqs. ([2]) and 

What is the origin of such an influx? Recent MD sim- 
ulations [ssl . [sgI . [stI of gas dissolved in water in con- 
tact with a surface demonstrated that on hydrophobic 
surfaces there is gas enrichment near solid walls. Such 
an enhanced gas concentration at hydrophobic walls has 
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n|; also neutron 



been confirmed spectroscopically 
ray reflectivity measurements [1^ 

reflectivity experiments [sll suggest it. In the MD simu- 
lations of ref. [3^ , for sufficiently hydrophobic walls the 
concentration near walls can exceed more than two or- 
ders of magnitude above the concentration in bulk liq- 
uid. Physically, this enrichment occurs because there is 
a potential (j){z) attracting solute molecules to the wall. 
The equilibrium concentration of the solute is thus de- 
termined by the balance between diffusion and attrac- 
tion, according to = — where C is the mo- 
bility of the solute. The Einstein relation D — ksT/^ 
then implies the equilibrium distribution of solute as 
c(z) = Co exp{—(j){z) / ksT) . A hundredfold concentration 
of solute molecules near the wall as found for the case 
studied in ref. [1^1 implies that (j){z) decreases by about 
AJ^bT near the wall, over a molecularly determined length 
scale. In general, we assume that the energy gain A0 of 
a solute molecule at the wall equals A(/) = sksT. 

Solute molecules at a gas-liquid interface are attracted 
to the wall and hence driven inside the bubble, see again 
figure 1. The size of the gas influx Ji„ is given by 



J„(i?) w 27r 



^ l#(z 
r dr c{r) 



hir)) 



C 



dz 



(4) 



Here we have again approximated the bubble shape h{r) 
as nearly flat. Since the attractive force drops off quickly 
with distance, the flux is only appreciable near the con- 
tact line. Using the fact that near the contact line 




FIG. 1: Sketch of gas outflux and influx into a surface 
nanobubble and definition of tlie contact angle 9, the surface 
nanobubble radius R, and the radius of curvature Re- 
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FIG. 2: Gas outflux jout (black-solid, from eq. ((Sjl) and influx 
jin (red-dashed, from eq. ((SJ) into the surface nanobubble as 
function of bubble radius R. The crossing point deflnes the 
equilibrium radius R* . If the slope of jout at ii = is larger 
than that of jin , no surface nanobubbles can emerge. For this 
illustrative plot we used eq. ((7)1 with 5 = 70nm, 9oo = 50°, 
and 9o = 0° , and the values s = 0.36 for the relative strength 
of the attraction potential, Coo/co = 0.25 for the relative gas 
concentration, and D = 10~^m^/s (which only scales the y- 
axis). For these data the stable equilibrium radius is R* = 
85nm. 



h(r) — [R — r) taii9 we can approximate equation ([¥]) to 
obtain the mass influx Jin{R) ~ 27ri?c(i?)A(/!)/(Ctan6'), 
or the volume influx 



jin (R) 



J„(i?) 2TrsDR 



c{R) tane{R)' 



(5) 



As anticipated above, eq. (O suggests that the in- 
flux scales linearly with R, exactly the same scaling 
as the diffusive outflux. However, recent measurements 
38 , 4C , Hll, \^ , including some on surface nanobubbles 
4J|, have shown that for sufficiently small bubbles 



43 



3 



150 



100 



DC 




FIG. 3: Dependence of equilibrium radius R* on the wall 
attraction strength s for relative gas concentrations Coo/co ~ 
0.1, 0.25, 0.5, 0.95, bottom to top. 



or droplets, the contact angle 9 depends on the contact 
line curvature and thus on the bubble size. This can be 
described through a modified Young-Dupre equation 



cos(7r — 6) ~ cos(7r — 6*00) 



C 



R + S' 



(6) 



where C and 6 are constants with the dimensions of 
length, and 6rx, is the usual macroscopic equilibrium con- 
tact angle. A popular (though disputed j39|) model as- 
sumes that the correction arises from a line tension t 
of the contact line [1^, in which case C = T/7, where 
7 is the liquid surface tension. The ratio C/S is set 
by the contact angle f?o for R ^ 0, namely C/S = 
cos(7r — 9oc) — cos(7r — Oq), so that one can rewrite ([6|) 
as 



cos 9 — cos 9^ 



cos 9r, 



cos 



1 + R/S 



(7) 



The length scale S is the onset-scale of microscopic correc- 
tions to the macroscopic contact angle. The exact form 
of eq. ((7]) and the exact values of the parameters are not 
relevant in the context of this paper; all what is needed 
is that 9{R) decreases with decreasing bubble radius R, 
which indeed is observed in refs. [s^ [40, [4^ [43] . We take 
S = 70nm, 9^ = 50°, and 6*0 = 0°. Note that the behav- 
ior 9{R ^ 0) = leads to a stabilization of small surface 
nanobubbles, as their curvature becomes small. In the 
limiting case i? ^ it even vanishes and the numerically 
found [35| case of a mono- or bilayer of gas molecules 
between the surface and the liquid is recovered. 

The outflux jout{R) and the influx jin{R) are shown 
in figure 2 for some typical parameters. In that case, 
a stable dynamic equilibrium radius R* defined through 
jout{R*) = jin{R*) exists. Note that indeed various stud- 
ies have revealed the existence of a prefered radius of the 
surface nanobubbles [1, |T^, which depends on the gas 
concentration p^ . For bubbles smaller than the equi- 
librium radius, the influx overcompensates the outflux. 



for larger bubble the outflux wins. If the gas concen- 
tration Coo decreases, the outflux increases and the equi- 
librium radius is becoming smaller, in agreement with 



13. If, 



on 



experimental observations on nanobubbles 
the other hand, the surface gets more hydrophobic and 
^oo decrease and/or the attraction potential s increases, 
the equilibrium radius is shifted towards larger values, 
again in agreement with experimental observations. 

The necessary condition for a stable dynamic equilib- 
rium and therefore for stable nanobubbles to exist is that 
at small bubble sizes R ^ the influx jin{R) is larger 
than the outflux jout{R)- This implies the condition 



s > - tan 9o 



c{R 0) 



— tan 9n 
2 



(8) 



for small gas concentration Coo ^ c{R — > 0). The condi- 
tion ([5]) is satisfied for sufficiently hydrophobic surfaces, 
where tan{9{R — > 0)) is small enough and when the solute 
is attracted to the wall, so that s is large enough. If con- 
dition is fulfilled, there also is an unstable equilibrium 
at _R* = 0, allowing for spontaneous nanobubble forma- 
tion. If that condition is not fulfilled, the equilibrium at 
R* = becomes stable and no surface nanobubbles form, 
in spite of an enhanced gas concentration very close to 
the surface. These may have been the conditions of the 
X-ray reflectivity experiments of ref. [Toj . where a gas 
accumulation, but no nanobubbles, have been seen. 

How does the prefered bubble radius R* depend on 
the material properties and the control parameters? To 
illustrate this dependence, we show i?*(s, Cqo/cq) result- 
ing from our simple model (eqs. ([31 and d?]), see fig. 
131 The dependence of R* on the attraction strength s 
is as expected rather strong - slight (chemical or struc- 
tural) inhomogeneities in the surface will therefore result 
into some distribution in the prefered radius, as indeed 
experimentally seen in refs. p, [l5| . In contrast, away 
from saturation, the dependence on the exact value of 
the relative gas concentration Coo/cq is relatively weak. 

Some comments on the driving mechanism: While the 
gas diffuses in the liquid, there is a gas flow from the con- 
tact line towards the surface inside the bubble. To avoid 
conflict with the second law of thermodynamics, such a 
state can only be transient. We do not know the origin 
of the non-equilibrium; it could be caused by tempera- 
ture gradients, chemical effects, or local oversaturation. 
In principle the observation method itself (e.g. AFM or 
optical detection of the nanobubbles) could be the origin 
of the equilibrium distortion. If the system were iso- 
lated, the second law of thermodynamics requires that 
in the long term (possibly hours or even days) the driv- 
ing mechanism would expire and the surface nanobubbles 
therefore dissolve, i.e., nanobubbles would be a transient 
phenomenon. Unfortunately, long-term observations of 
surface nanobubbles in closed sytems have to our knowl- 
edge not yet been carried out. Note that the required 
non-equilibrium situation must also reflect in the chemi- 
cal potential ^, as in equilibrium gas flow can only occur 



4 



along a gradient in fj,, and a circular flow would not be 
possible. 

A method to create controlled non-equilibrium condi- 
tions is to generate surface nanobubbles through elec- 
trolysis [111, H^. In ref. ^461] it has been shown that 
several tens of seconds after switching on the potential, 
these surface nanobubbles do not further grow, in spite 
of a nonzero current. This observation suggests that the 
nanobubbles indeed have achieved a dynamic stable equi- 
librium, with the Laplace pressure driven gas outflux be- 
ing compensated by the gas influx at the electrode. 

In summary, we have suggested a dynamic equilibrium 
stabilization mechanism for surface nanobubbles: The 
gas outflux driven by Laplace pressure is compensated by 



a gas influx at the contact line, leading to a metastable 
equilibrium under certain conditions. A necessary ingre- 
dient into this model is a contact angle decrease with 
decreasing bubble size R. Though such a decrease has 
been observed both experimetally [H, and numeri- 
cally [H, SO] , further work is required to better quantify 
and understand this dependence d{R). 
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